Research in Astron. & Astrophys. Vol.0 (200x) No.O, 000-000 „ , . 

Research in 

(http://www.raa-journal.org) Astronomy and 

Astrophysics 



Metal Abundance and Kinematical Properties of M81 Globular 
Cluster System 

Jun Ma^'^ *, Zhenyu Wu^ Tianmeng Zhang^ Song Wang^'^, Zhou Fan', Jianghua Wu^, Hu 
Zou', Cuihua Du^, Xu Zhou' and Qirong Yuan'' 

'National Astronomical Observatories, Chinese Academy of Sciences, Beijing, 100012, P. R. China 
^Key Laboratory of Optical Astronomy, National Astronomical Observatories, Chinese Academy of 
Sciences, Beijing, 100012, China 

^Graduate University, Chinese Academy of Sciences, Beijing, 100039, China 
''Department of Astronomy, Beijing Normal University, Beijing 100875, China 

^College of Physical Sciences, Graduate University of the Chinese Academy of Sciences, Beijing 100049, 
China 



Department of Physics, Nanjing Normal University, Wen Yuan Road 1, Nanjing 210046, China 



Abstract In this paper, we presented metal abundance properties of 144 M81 globular clus- 
ters. These globulars consist of the largest globular cluster sample in M81 till now. Our 
main results are: the distribution of metallicities are bimodal, with metallicity peaks at 
[Fe/H] «i -1.51 and -0.58, and the metal-poor globular clusters tend to be less spatially 
concentrated than the metal-rich ones; the metal-rich globular clusters in M81 do not demon- 
strate a centrally concentrated spatial distribution as the metal-rich ones in M31 do; like our 
Galaxy and M31, the globular clusters in M81 have a small radial metallicity gradient. These 
results are consistent with those obtained based on a small sample of M8 1 globular clusters. 
In addition, this paper showed that there is evidence that a strong rotation of the M8 1 globu- 
lar cluster system around the minor axis exists, and that rotation is present in the metal-rich 
globular cluster subsample, while the metal-poor globular cluster subsample shows no evi- 
dence for rotation. The most significant difference between the rotation of the metal-rich and 
metal-poor globular clusters occurs at intermediate projected galactocentric radii. The results 
of this paper confirm the conclusion of Schroder et al. that M8 1 's metal-rich globular clusters 
at intermediate projected radii were associated with a thick disk of M8 1 . 
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1 INTRODUCTION 

An understanding of galaxy formation and evolution is one of the principal goals of modern astrophysics. 
As the oldest stellar systems in the universe, globular clusters (GCs) maintain a fossil record of the early 
history of galaxies, so they are important tracers of galaxy formation and evolution. In addition, since GCs 
are very dense, gravitationally bound spherical systems of several thousands to more than a million stars, 
they can be observed out to much greater distances than individual stars. 

The metallicity (or color) distribution of GCs is of particular importance in deepening our knowl- 
edge of the dynamical and chemical evolutio n of the parent galaxies. For example, based o n data from 
the Hubble Space T elescope (HST) archive, iGebhardt & Kissler-Patigl (Il999h . iLarsen et alJ (l200lh and 
iKundu & Whitmord (l200lb showed that many large galaxies possess two or more subpopulations of GCs 
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that have quite different chemical compositions (see also IWest et all I200I . IPeng et alJ (l2006h presented 
the color distributions of GC systems for 100 early-type galaxies from the Advanced Camera for Surveys 
(ACS) Virgo Cluster Survey, and found that, on average, galaxies at all luminosities appear to have bimodal 
or asymmetric GC color/metallicity distributions. The presence of color bimodality among these old GCs 
indicates that there have been at least two major star-forming mechani sms in the (early) histories of massive 
galaxies (IWest et al.ll2004lPeng et al.ir200alBrodie & Stradeiil2006l) . 

' Fan et alj (|2008') presented a metallicity distribution of 148 Galactic GCs (iHarrisl [l996l 2003 updated 
version) that apparently shows two peaks (i.e., two distinct metal-poor and metal-rich GC populations). A 
double-Gaussian can best fit these two subpopulations, t he mean metallic ity values are - 1 .620 and -0.608 
dex, respectively. Using the data for 431 GCs in M31. 'L ee et al.l (l2008l) studied the metallicity distribu- 
tion, which is asymmetric, implying the possibility of bimodality. The bimodal test results in the two 
group : the metal-poor \yith a mean value of -1.47 and the r netal-rich with a mean value of -0.62 (see 
also Barmbv et aUl2000HPe"rrett et am2002l:lFan et an.l2008l) . lMa et alj (l2005h presented a metalhcity dis- 
tribution of 94 M81 GCs being bimodal, with metallicity peaks at -1.45 and -0.53. 

M81 is one of the nearest Sa/Sb-type spiral galaxies outside the Local Group, very similar to M31, and 
roughly as massive as the Milky Way. So, beyond the Local Group, it is a good candidate for reaching a 
detailed study of a spiral gala xy GC system for compa rison to the Milky Way and M3 1 GC system. In the 
pioneering work of M8 1 GCs. lBrodie & Huchral ( ll99T h derived spectrosc opic metallicities for eig ht GCs in 
M8 1 and presented the mean of [Fe/H] = -1.46 + 0.31 for them. Then, Perelmuter et alj d 19951) obtained 
low signal-to-noise spectra of 82 GC candidates, 25 of which were confirmed as bona fide M81 GCs. 
These authors derived the mean metallicity to be [Fe/H] — —1.48 + 0.19 both from the weighted mean of 
the individual metallicities, and directly from the composite spectrum of these 25 GCs. To maximize the 
success rate of M8 1 GC candidate list for the ongoing spectroscopic observations, i Perelmuter & Racing 
used an extensive database that included photometric, astrometric, and morphological information 
on 3774 objects covering over a > 50' diameter field centered on M81 to reveal 70 GC candidates. 

ISchroder et al.' ('2002) presented moderate-resolution spectroscopy for 16 GC candidates from the GC 
list in Perelmuter & Racine (1995), and confirmed the se 16 candida tes as bona fide GCs. They derived 
metallicities for 15 of the 16 GCs. From their results. ISchroder et all (I2OO2I) concluded that the M81 GC 
system is very similar to the Milky Way and M31 systems, both chemically and kinematic ally. 

With the superior resolution of the HST, M8 1 is close enough for its clusters to be easily resolved 
on the basis of image structure (Chandaret al., 2001). Thus, using the B, V, and / bands of HST Wide 
Field Planetary Camera 2 (WFPC2), Chandar et al.1 (12001 ) imaged eight fields covering a total area of 
~ 40 arcmin^, and detec ted 1 14 compact star clusters in M81, the 59 of which are GCs. Based on the 
estimated intrinsic colors, IChandar et al.l (|2004|) found that the M8 1 GC system has an extended metallicity 
distribution, which argues the presence of both metal-rich and metal-poor GCs. 

iMa et all ( [20051 120061 l2007h presented detailed studies on M81 GCs based on the spectral and multi- 
color observations such as: the distributions of intrinsic B and V colors and metallicities of 95 GCs; the 
spectral energy distributions of 42 GCs in 13 intermediate-band filters from 4000 to 10,000 A, using the 
CCD images of M81 observed as part of the Beijing- Arizona-Taiwan-Connecticut (BATC) multicolor sur- 
vey of the sky; the spatial and metal abundance properties of 95 GCs. 

iNantais et al.l ( l2010bl) presented a catalog of exte nded objects in and around M8 1 based on an HST ACS 
/-band mosaic, and found 233 good GC candidates. iNantais & Huchral (l2010l) obtained spectra of 74 GCs 
in M8 1 wit h Hectospec on the 6.5 m MMT on Mt. Hopkins in Arizona, and determined their metallicities. 
INantais &~H uchra (2010) analyzed the kinematics and metallicity of the M81 GC system combined their 
results with 34 M81 GC velocities and 33 metallicities from the literature. 

In this paper, we will analyze the metallicity of the M8 1 GC system based on a large sample of GCs, 
and analyze the kinematics of the M81 GC system in details. The outline of the paper is as follows. In 
§ 2, we provide the analysis of the metallicity, and we provide the analysis of the kinematics in § 3. The 
summary is presented in § 4. 
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Fig. 1 Positions of th e 144M81 GCs. The large ellipse is the D25 boundary of the M81 disk from 
Ide Vaucouleurs et alj ( Il99lh . 



2 METALLICITY ANALYSIS 



2.1 A sample of GCs 



iNantais & Huchral (1201 Ol) obtained spectra of 74 GCs in M81 with Hectospec on the 6. 5 m MMT on M t. 
Hopki ns in Arizona, and determined their metallicities using a method based on the iBrodie & Huchral 
([l9W) presc ription, but w ith new index versus metallicity calibrations determined from the 41 Milky Way 
spectra from lSchiavon et al. (2005), degraded to match the 5 A resolution of t heir spectra (see J Mantais et all 
|2010a, for details). In order to increase the sample size and spatial extent. iNantais & Huchral (120101) in - 
clu ded the derived metalhcit ies of 33 GCs in M81 from I Schroder et alj ( 12002'). 'Perel muter et alj ( Il995h . 
and iBrodie & Huchral ( Il99lh . Altogether, the re are 107 GCs in the INa n tais & Huchral ( 1201 Ol) sample of 
M81 GCs for performing metallicity analysis. Ma et al. (2005) studied the distributions of metallicities of 
95 M81 GCs which are from .Perelmuter et al. ( 1995). Chandar et al. (.2001.) and Schroder et al. C2002.) . By 
cross-checking the sample GCs of between |Ma et a il (I2005I) and INantais & Huchral (l2010l) . there are 51 
GCs in common. We adopt'Nantais & Huchra' ('2010') metallicities with uncertainties smaller than 1.0 dex 
whenever available. However, when Nantais & Huchra (2010) metallicity uncertainties are la rger than 1.0 
dex, we adopt the metallicities collected b y iMa et al] (l2005l) . In addition, as iMa et aP (l2005h pointed out 
that, since GC 96 of IChandar eTal] (12001*) has very high {B - V)o {{B - V)o = 1.778), and the metallicity 
obtained using the color-metalhcity correlation is too rich (0.95 dex), we do not include it in our sample of 
M81 GCs. Altogether, there are 144 M81 GCs in our sample. 
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Fig. 2 Histogram of the metallicities of 144 M81 GCs. 

Figure 1 shows the positions of 144 GCs in M81 with respect to the M81 disk. The X coordinate is the 
position along the major axis of M81, where positive X is in the northeastern direction, while the Y coordi- 
nate is along the minor axis of the M81 disk, increasing towards the northwest. The relative coordinates of 
the M81 GCs are derived b y assuming standard geometric pararneters f or M81. The distance modulus for 
M81 is adopted to be 27.8 (iFreedman et al.lll994l; IChandar et all l200l'). We adopted a central position for 
M81 atao = 09*'55™33^780 and (5o = +69°03'54".34 (J2000.0) following , Schroder et al.. (,2002). Formallv. 

X^Asm6 + Bcos6 ; (1) 

y = -Acos6l + Bsin6i , (2) 

where A = sin(Q' - q-q) cos 6 and B = sin 6 cos dp - cos(a - ao) cos 6 sin Sq. We adopt a position angle of 
= 157° for the major axis of M81 dChandar et al.Ll200lh . 

2.2 Metallicity distribution 

Figure 2 displays the metallicity of 144 GCs in M81, and two metallicity peaks appear clearly. To ma ke 
quantitative statements about the bimodal metallicity distribution, a KMM test dAshman et al.l 1 19941) is 
applied to the data. This test uses a maximum likelihood method to estimate the probability that the data 
distribution is better modeled as a sum of two Gaussians than as a single Gaussian. Here we use a ho- 
moscedastic test (i.e., the two Gaussians are assumed to have the same dispersion). The metallicities of the 
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Fig. 3 Spatial distributions of the metal-rich and metal-poor GCs in M8 1 . 



two peaks, the f -value, and the numbers of GCs assigned to each peak by the KMM test are [Fe/H] -1.51 
and -0.58, 0.02, and 91 and 53. The f-value is in fact the probability t hat the data are d rawn from a single 
Ga ussian distribution. It is evi dent that our results are in agreement with lMa et al.l (l2005h (- 1 .45 and -0.53) 
and lNantais & Huchral(l2010l) (-1.55 and -0.61). 

2.3 Spatial distribution 

Figure 3 shows the projected spatial distributions of the metal-poor and metal-rich GCs in M81. From 
Figure 3, it is cl ear that the metal-r i ch GCs in M81 are not as centrall y concentrated as the metal-rich 
GCs of M31 are (iHuchraet al.lll99ll;|Pen-ett et al.ll2002l:lFan et al.Ll2008l) . and the metal-poor G Cs tend to 
be les s spatially concentrated than the metal-poor GCs. These results confirm the conclusion of i Ma et al] 
(l2007h based on the small sample of M81 GCs. 

Figure 4 presents the histogram for the metal-poor and metal-rich GCs in M8 1 . It shows that most of the 
metal-rich GCs distribute at projected radii of 1-7 kpc, and the peak lies at 2-4 kpc. It is also noted that the 
metal-rich GCs distribute within the inner 13 kpc, while the metal -poor GCs do out to radii of ~ 20 kpc. In 
the Milky Way, the metal-rich GCs reveal significan t rotations and have hist orically been associated with the 
thick-disk system (IZinnlll985l; lArmandroffl Il989l) : however, other works jFrenk & Whitel 1 1 982t iMinnitil 



ll995l:IC6tel[l999[lForbes et al.ll200lb suggested that metal-rich GCs within ~ 5 kpc of the Galactic center 
are better associated with the bulge and bar. In M31, Elson & Walterbos (1988) showed that the metal- 
rich GCs c onstitute a more hig hly flatted system than the metal-poor GCs, and appear to have disklike 
kinematics:'Huchra et al.l d 199 lb s howed that the metal-rich GCs are preferentially close to the galaxy center. 
In addition, Huchra et al] (1 199 lb showed that the distincti on between th e rotation of the metal-rich and 
metal-poor GCs is most apparent in the inner 2 kpc. So, iHuchra et al.l (il99lb concluded that the rich- 
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Fig. 4 Radial distributions of the metal-rich and metal-poor GCs in M8 1 . 



metal GCs in M31 appear to form a central rotating disk system. With a large sample of 321 velocities, 
[Perrett et al. (2002) provided a comprehensive investigation on the kinematics of the M3 1 GC system. These 
authors showed that the metal-rich GCs of M3 1 appear to constitute a distinct kinematic subsystem that 
demonstrates a centrally spatial distribution with a high rotation amplitude, but do not appear significantly 
flattened, which is consistent with a bulge population. It is of interests to me ntion that , Schroder et al. ( 2002) 
performed a maximum-likelihood kinematic analysis on 166 M31 GCs of iBarmby et all (120001) and found 
that the most significant difference between the rotati on of the metal-rich an d metal-poor GCs occurs at 
intermediate projected galactocentric radii. Especially. ISchroder et al. (12002) pres ented a potential thick- 
disk population among M31's metal-rich GCs. For M81 GCs, Schro der et al. ( 2002*) performed a kinematic 
analysis of the velocities of 44 M81 GCs, and strongly suggested that the metal-rich GCs are rotating in the 
same sense as the gas in the disk of MS 1 . ISchroder et aU (120021) concluded that, although their GC sample 
is not large enough to make a direct comparison between metal-rich and metal-poor GCs in specific radius 
ranges, the conclusion with MS l's metal-rich G Cs at intermediate projected radii being associated with a 
thick disk in M81 is correct. In ISchroder et al.l (|2002), their sample is not large enough to make a direct 
comparison between metal-rich and metal-poor GCs in specific radius range. In this pape r, with a larger 
GC sample, we can do this comparison (see §3). Our results will confirm the conclusion of ISchroder et al.l 
(l2002h that MSl's metal-rich GCs at intermediate projected radii is associated with a thick disk in M81. 



Globular Clusters in M81 



7 



-| , , , , 1 , , , , 1 , , , , 1 , , , , 1 . , r- 

Full Sample 




0) 





J u_ 



Metal-Poor 



Metal-Rich 



5 10 15 20 

Projected Radius (kpc) 



25 



Fig. 5 Metallicity as a function of projected radius for M81 GCs. 
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Fig. 6 Positions of th e 108M81 GCs. The large ellipse is the D25 boundary of the M81 disk from 
Ide Vaucouleurs et alj ( Il99lh . 



2.4 Metallicity gradient 

One early formation model bv lEggen et alj (Il962h argued for a single, large-scale collapse of material to 
form galactic bodies such as the Milky Way, in which the enrichment timescale is shorte r than the col- 
lapse time, the halo stars and GCs should show large-scale metallicity gradients; however. ISearle & ZinnI 
(Il978i) presented a chaotic scheme in the early evolution of a galaxy, in which loosely bound pre-enriched 
fragments merge with the main body of the proto-galaxy over a significant period, so there should be ho- 
mogeneous metallicity di stribution. 

For the Milky Way, lArmandro3 (Il989h showed some evidence that metallicity gradients with both 
distance from the Galactic plane and distance from the Gala ctic center were present in the disk GC system. 
For M31, there are some inconsistent conclusions, such as Ivan deii Berghl(ll969h who showed that there 
is little or no evidence for a correlation between metallicity a nd projected ra dius, but most of his GCs 
were inside 50"; however, some authors (see, e.g. iHuchra etaU .1982; Sharovi Il988l; iHuchra et all Il99ll 
iBarmbv et al.ll2006l;|Perrett et al.ll2002h pres e nted e vidence for a weak but measurable metallicity gradient 
as a function of projected radius. lpan et al. I (I2OO8') confirmed the latter conclusion based on their large 
sample of spectral metallicity and color-derived metallicity. Figure 5 plots the metallicity of the M8 1 GCs 
as a function of galactocentric radius. Clearly, the dominant feature of this diagram is the scatter in metallity 
at any radius. Since our sample GCs are mainly distributed in the inner 10 kpc, it is difficult to determine the 
metallicity gradient. It is true that, smooth, pressure-supported collapse models of galaxies are unlikely to 
produce a result like this. In order to present a quantitative conclusion, we made least-squares fits: the total 
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Fig. 7 Errors of the measured velocities of M8 1 GCs versus a galactocentric radial distance. 



sample o f GCs does not have a significant metallicity gradient (-0.016 + 0.012 dex kpc which confirm 
Ma et al. (2007) conclusion based on a smaller sample of M81 GCs. This result is also in agreement with 
Kong et al., C2000.) . who o btained metallicity maps o f M8 1 field by comparing simple stellar population 
synthesis models of BC96 (iBruzual & Charlol 1 19961) with the integrated photometric measurements of the 
BATC photometric system, and did not find, within their errors, any obvious metallicity gradient from the 
central region to the bulge and disk of MS 1 . In addition, we should emphasize that, in the least-squares fits 
of this paper, the metal-rich GCs seem to act an important part in determining the metallicity gradient. 

3 KINEMATICS 

The M81 GC sample for studying kinematics is from iNantais & Huchral ( 1201 Ol) . who obtained spectra of 
a total of 207 extended objects in M81 with Hectospec on the 6.5 m MMT on Mt. Hopkins in Arizona. 
These objects were selected from iNantais et al] (l2010 bV who presented a catalog of extended objects in 
the vicinity of M81 based on a set of 24 HST AC S Wide Field Came r a (WF C) F814W images, and found 
233 good GC candidates. Based on these spectra. INantais & Huchral (I2OIOI) fou nd 74 GCs in M8 1 , 62 o f 
which were newly confirmed by themselves. Combined with other 34 GCs from B rodie & Huchral (Il99lb . 
|Perelmuter et al. ( 1995), and Schroder et al. (2002), Nantai s & Huchraj 2 C)10) p resented a catalog of radial 
velocities for 108 GCs in M81. Based on this catalog. INantais & Huchral (l2010l) studied the rotation, mean 
velocity, and velocity dispersion of M8 1 GCs, and presented a conclusion that, the M8 1 GC system as a 
whole shows strong evidence of rotation, with V,- = 108 + 22km s"' (a deprojected rotational velocity) 
overall. In this paper, we will present other evidence of rotation of the M8 1 GC system. 
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Fig. 8 Radial velocity histogram for 108 GCs in M81. The vertical dashed lines indicates the 
mean value of the radial velocities, (v) = -24 km s ' . 



Figure 6 shows the positions of 108 GC s in M8 1 with respect to the M81 disk. Using the maximum 
likelihood method o f |Pryor&Meyianl(ll993h (see also lWu et alll2009l for details), we find a mean line-of- 
sight velocity for the full GC sample of (v) - -24 +14 km s"\ somewhat higher than the systemic velocity 
of M31, (v) - -34 + 4km s"' (de Vaucouleurs et al. 1991). The GCs have an overall velocity dispersion 
of o",, - 142 +11 km s"'. In addition, the mean value of the radial velocities and the velocity dispersion 
are derived to be ( v) = -22 ± 25 km s"' and c^^, = 153;^[g km s"' using the biweight location and scale of 
iBeers etalj(ll990l) . The kinematics of metal-rich and metal-poor GC populations will be discussed later in 
this section. First we examine the global kinematic properties of the M81 GCs. Figure 7 shows the errors of 
the measured velocities of M8 1 GCs versus a galactocentric radial distance. From Figure 7, we cannot find 
evidence that the mean errors for the measured velocities are larger in the inner region where background 
level is higher than that in the outer region of M81. Figure 8 displays the radial velocity histogram for 108 
GCs. It appears that the velocity distribution for all GCs is roughly symmetric with respect to the mean 
value of the radial velocities 

Figure 9 presents the spatial distribution of M8 1 GCs with the measured velocity. Here we adopted as 
the sys tematic velocity of M3 1 the value derived using the maximum likelihood method of Prvor & MeylaiJ 
(Il993h . (v) = -24± 14 km s"' . In Figure 9, we present the radial velocities for M8 1 GCs versus the projected 
radii along the major axis (X) and minor axis (Y), respectively. A linear least-squares fit, passing through 
{X, V - {v))-{0, 0), along the major axis results in v - (v) - (-1 1.8 ± 2.7) X km s"', and passing through 
(Y, V - (v))=(0, 0), along the minor axis results in v - (v) = (-0.38 ± 5.0) Y km s"', where X and Y are units 
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Fig. 9 Radial velocity of the M8 1 GCs versus projected radius along the major axis {left) and the 
minor axis {right). Solid lines indicate the linear least-squares fits. 



of kpc. The large value of the slope between v - (v) and X indicates a strong rotation of the M81 GC system 
around the minor axis, while the small value of the slope between v - (v) and Y indicates no significant 
rotation around the major axis. 

Below, we will look for evidence for differences in subpopulations in M81 GCs that differ in their 
spatial and chemical characteristics. The formula for the rotation curve is v = Vsys + vvot sin(0 + 6*0), where 
v'sys is the mean velocity of the GC system. The results of our kinematics analysis are summarized in Table 
1. Figure 10 illustrates that rotation is present in the metal-rich GC subsample, while the metal-poor GC 
subsample shows no evidence for rotation. Table 1 also shows that the most significant difference between 
the rotation of the metal-rich and metal-poor GCs occurs at intermediate projected galactocentric radii. In 
addition, most of the metal-poor GC subsample have poorly constrained rotation axes and have rotation 
velocity-to-velocity dispersion ratios that are consistent with their being supporte d primarily by therma l 
motion rather than by rotation. These results are consistent with the conclusions of ISchroder et al.l (|2002|) 
forM31 GCs. 

4 SUMMARY 

In this paper, we presented metal abundance properties of 144 M81 GCs which consist of the largest GC 
sample in M81. We also studied the kinematics properties of the M81 GC system. Our main conclusions 
are as follows: 

1. The distribution of metallicities are bimodal, with metallicity peaks at [Fe/H] « -1.51 and -0.58. 

2. The metal-poor GCs tend to be less spatially concentrated than the metal-rich ones. 




Fig. 10 Quantitative assessment of rotation in M8 1 metal-rich (MR) and metal-poor (MP) GC 
subsamples based on projected galactocentric radius (Rp). 



3. The metal-rich GCs in M81 do not demonstrate a centrally concentrated spatial distribution as the 
M31 metal-rich GCs do. 

4. The GCs in M8 1 have a small radial metallicity gradient. 

5. There is evidence that a strong rotation of the M81 GC system around the minor axis exists. 

6. Rotation is present in the metal-rich GC subsample, while the metal-poor GC subsample shows no 
evidence for rotation, and the most significant difference between the rotation of the metal-rich and metal- 
poor GCs occurs at intermediate projected galactocentric radii. 

7. M8rs metal-rich GCs at intermediate projected radii are associated with a thick disk in M81. This 
result confirms the conclusion o f ISchroder et ai](l2002l) . 

ACKNOWLEDGMENTS 

We would like to thank the anonymous referee for providing rapid and thoughtful report that helped im- 
prove the original manuscript greatly. This work was supported by the Chinese National Natural Science 
Foundation grands No. 10873016, and 10633020, 10803007, 11003021, 11173016 and 11073032, and by 
National Basic Research Program of China (973 Program), No. 2007CB815403. 

References 



Armandroff, T E. 1989, AJ, 97, 375 



Globular Clusters in M81 



13 



Ashman, K. A., Bird, C. M., & Zepf, S. E. 1994, AJ, 108, 2348 

Barmby, P., Huchra, J., Brodie, J., Forbes, D., Schroder, L., & Grilhnair, C. 2000, AJ, 119, 727 

Beers, T. C, Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 

Brodie, J. P., & Huchra, J. P 1990, ApJ, 362, 503 

Brodie, J. P, & Huchra, J. P 1991, ApJ, 379, 157 

Brodie, J., & Strader, J. 2006, ARA&A, 44, 193 

Bruzual, G., & Chariot, S. 1996, unpubhshed 

Chandar, R., Ford, H. C, & Tsvetanov, Z. 2001, AJ, 122, 1330 

Chandar, R., Whitmore, B., & Lee, M. G. 2004, ApJ, 61 1, 220 

C6te,R 1999, AJ, 118, 406 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H., Buta, R., Paturel, G., & Fouque, P. 1991, Third 

Reference Catalog of Bright Galaxies, Vol. 2 (New York: Springer), 548 
Eggen, O. J., Lynden-Bell, D., & Sandage, A. R. 1962, ApJ, 136, 748 
Elson, R. A., & Walterbos, R. A. M. 1988, ApJ, 333, 594 
Fan, Z., Ma, J., de Grijs, R., & Zhou, X. 2008, MNRAS, 385, 1973 
Forbes, D. A., Brodie, J. P, & Larsen, S. S. 2001, ApJ, 556, L83 
Freedman, W. L., Wilson, C. D., & Madore, B. F 1994, ApJ, 427, 628 
Frenk, C. S., & White, S. D. M. 1982, MNRAS, 198, 173 
Gebhardt, K., & Kissler-Patig, M. 1999, AJ, 118, 1526 
Harris, W. H. 1996, AJ, 112, 1487 

Harris, W. E., Whitmore, B. C, Karakla, D., Okon, W, Baum, W A., Hanes, D. A., & Kavelaars, J. J. 2006, 
ApJ, 636, 90 

Huchra, J. P, Brodie, J. P, & Kent, S. M. 1991, ApJ, 370, 495 
Huchra, J., Stauffer, J., & van Speybroeck, L. 1982, ApJ, 259, L57 

Kong, X., et al. 2000, AJ, 1 19, 2745 

Kundu, A., & Whitmore, B. C. 2001, AJ, 121, 2950 

Larsen, S. S., et al. 2001, AJ, 121, 2974 

Lee, M. G., Hwang, H. S., Kim, S. C, Park, H. S., Geisler, D., Sarajedini, A., & Harris, W. E. 2008, ApJ, 

674, 886 
Ma, J., et al. 2005, PASP, 1 17, 256 
Ma, J., et al. 2006, PASP, 1 18, 98 
Ma, J., et al. 2007, PASP, 119, 1085 
Minniti, D. 1995, AJ, 109, 1663 
Nantais, J. B., & Huchra, J. P 2010, AJ, 139, 2620 

Nantais, J. B., Huchra, J. P, Barmby, P, & Olsen, K. A. 2010a, AJ, 139, 1178 

Nantais, J. B., Huchra, J. P, McLeod, B., Strader, J., & Brodie, J. P 2010b, AJ, 139, 1413 

Peng, E. W., et al. 2006, ApJ, 639, 95 

Perelmuter, J. M., Brodie, J. P, & Huchra, J. 1995, AJ, 110, 620 
Perehnuter, J. M., & Racine, R. 1995, AJ, 109, 1055 
Perrett, K. M., et al. 2002, AJ, 123, 2490 

Pryor, C, & Meylan, G. 1993, in ASP Conf. Ser. 50, Structure and Dynamics of Globular Clusters, ed. S. 

G. Djorgovski & G. Meylan (San Francisco: ASP), 357 
Schiavon, R. P, Rose, J. A., Courteau, S., & MacArthur, L. A. 2005, ApJS, 160, 163 
Schroder, L. L., Brodie, J. P, Kissler-Patig, M., Huchra, J. P, & PhilUps, A. C. 2002, AJ, 123, 2473 
Searle, L., & Zinn, R. 1978, ApJ, 225, 357 
Sharov, A. S. 1988, Soviet Astron. Lett., 14, 339 
van den Bergh, S. 1969, ApJS, 19, 145 

West, M. J., Cote, P, Marzke, R. O., & Jordan, A. 2004, Nature, 427, 31 
Wu, Z. Y., Du, C. H„ Ma, J., & Zhou, X. 2009, ChPhL, 26, 9701 
Zinn, R. 1985, ApJ, 293, 424 



JunMa 

Table 1 Kinematics of M81 globular clusters. 



Spatial Group 


(km/s) 




(km/s) 




Vrot 

(km/s) 


00 
(deg) 


number 


Rp < 4kpc 
















MR 


-50.85 ± 31.76 


152, 


.07 ± 22, 


.80 


128 ± 27 


102.50 ± 19.05 


24 


MP 


-5.85 ±30.19 


149, 


.64 + 21 


.51 


40 ±56 


0.92 ± 69.90 


26 


2kpc < Rp < 6kpc 
















MR 


-39.02 + 28.72 


148, 


,32 + 28, 


.72 


141 + 15 


119.00 ± 18.16 


28 


MP 


-69.12 + 20.69 


111, 


,89 + 14, 


.95 


35 + 28 


-8.20 ± 54.86 


33 


3kpc <Rp< 7kpc 
















MR 


-31.77 ± 30.23 


138, 


.27 ± 22, 


.37 


138 ± 18 


118.60 ± 13.06 


22 


MP 


-29.52 ± 27.08 


127, 


.30 ± 20 


.81 


29 ±40 


184.78 ± 163.35 


25 


4kpc < Rp < 8kpc 
















MR 


-40.97 + 33.98 


127, 


.58 + 25, 


.49 


84 ±21 


116.01 ±43.54 


15 


MP 


-43.15 + 32.67 


133, 


,62 ± 25, 


.43 


21 ±39 


125.99 ± 340.47 


19 


Rp > 8kpc 
















MR 


-34.03 ± 48.97 


104, 


.72 ± 42, 


.02 


104 ±27 


119.92 ±51.11 


5 


MP 


6.06 ± 36.57 


125, 


.83 ± 28 


.42 


67 ±42 


92.59 ± 90.44 


14 



